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' With the Schwinger model as example I discuss properties of lattice Dirac oper- 

' ators, with some emphasis on Monte Carlo results for topological charge, chiral 

i fermions and eigenvalue spectra. 

1 Introduction 

\Q ' 1.1 Continuum Concepts vs. Lattice Concepts 

Relativistic particle physics is described by relativistic quantum field theory; 
£SJ ' the theories explaining the fundamental interactions are gauge theories. QFT 

as it is formulated in terms of Lagrangians and functional integration has to 
be regularized. The only known regularization scheme that retains the gauge 
symmetry is replacement of the space-time continuum by a space-time lattice. 
For convenience and other reasons this is done in an Euclidean world. 
| The most prominent QFT is QCD, the SU(3) gauge theory of quarks 

0^ ' and gluons. There various non-perturbative phenomena appear intertwined: 

confinement and chiral symmetry breaking. The classical continuum gauge 
fields A that are continuous and differentiable, living on compact manifolds, 
may be classified by a topological quantum number (the Pontryagin index) 
Q(A). Changing through continuous deformations of the field from one such 
| topologically defined sector to another is impossible. 

' Topology is closely related to the fermion zero modes; these are eigenstates 

(eigenvalue E = 0) of the Dirac operator 



><3 ■ i*fr(dr-ieAJtl> = Efl> = -> 75^ = ±V>- (1) 

H ' 

Due to the anti-commutation of 75 with T> the zero modes can be chosen 
as eigenstates of 75 with definite chirality. The Atiyah-Singer index theorem 
(ASITp relates the topological charge Q(A) of the background gauge field to 
these modes, 

Q(A) = index(A) = ri+ - n_ . (2) 
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where n± denotes the number of independent zero modes with positive or 
negative chirality, respectively. In two dimensions, another theorem of the 
continuum - the so-called Vanishing Theoremcl - ensures that only either 
positive or negative chirality zero modes occur. 

Quantization involves summation over non-differentiable fields and the 
lattice formulation does not provide for a unique definition of topological 
charge at all. Any lattice definition involves implicit or explicit assumptions, 
usually on the continuity and smoothness at the scale below one lattice spac- 
ing. Quantization usually means Monte Carlo integration on finite lattices 
and it is legitimate to question the ergodicity of the respective "simulation" 
with regard to the topological sectors. 

Even the task of putting fermions on the lattice introduces limitations. 
Various lattice Dirac operators have been proposed. It has been demonstrated 
earfytl that within quite general assumptions (like locality and reflection pos- 
itivity) it is not possible to have single chiral fermions; chiral symmetry has 
to be broken explicitly. For the simple Wilson action, the breaking of the 
symmetry is so bad that no trace of the chiral properties of the continuum 
theory is kept in the lattice theory. If, however, chiral symmetry is broken 
already on the level of the action, how can one hope to identify the (expected) 
spontaneous chiral symmetry breaking of the full theory? Is there a lattice 
version of the ASIT? As we will discuss below, recent developments do al- 
low us to construct lattice Dirac operators, which break chiral symmetry in a 
minimal way. 

Confinement, on the other hand, proved to be more straightforward. The 
gauge coupling g enters the lattice gauge action in form of the multiplicative 
coupling (3 = l/g 2 . The lattice formulation works particularly well at strong 
coupling (small 0), where confinement within a non- vanishing domain of the 
couplings was proved. In the Monte Carlo calculations for lattice QCD up 
to now no signal was found that indicates a phase transition between that 
confinement phase and the weak coupling (perturbative) regime; for finite 
temperature such a deconfinement transition was established. 

In order to retrieve continuum QFT numbers for the physical quantities 
- like masses of hadronic bound states or certain matrix elements - one has 
to show that all dimensional physical quantities scale according to the scaling 
function of the lattice spacing a(/3), which in turn should asymptotically (in 
the continuum limit f3 — > oo, a — > 0) agree with continuum renormalization 
group scaling. 

Current lattice studies have to live in that environment: Try to improve 
scaling properties in order to get reliable continuum results and try to deal 
with chirality without destroying it from the begin! 
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1.2 Lattice Dirac Operators 



All Dirac operators have 75-hermiticity, 

l5 v l5 = vK (3) 

It follows that 

• The eigenvalues are either real or are complex conjugate pairs. 

• The operator 75 T> has a real spectrum. 

• We denote by Vi the eigenvector of D for eigenvalue A^; then 75 Vi is an 
eigenvector of T> < for the same eigenvalue. 

• The diagonal entries of the chiral density matrix vanish for non-real eigen- 
values: Aj ^ Re — > («i|75 Vi) = 0. 

• The non-diagonal entries of the chiral density matrix vanish whenever 
their respective eigenvalues are not complex conjugate pairs: A; ^ Xj — » 
(Vi\j5Vj) = 0. 

An important conclusion is that only real eigenvalues lead to contributions 
to the diagonal elements. Thus these modes are the only candidates for zero 
modes in the continuum limit: in the continuum theory the (normalized) zero 
modes contribute to the diagonal entries, (^(75 fp) = ^l- 

Chirality is necessarily broken for lattice Dirac operators. However, barely 
noticed for almost two decades, Ginsparg and WilsorO formulated a condition 
(GWC) under which circumstances^ remnants of the chiral symmetry survive 
in a lattice action for massless fermions. If the Dirac operator obeys 

~ { l5 ,V} = aV l5 RV , (4) 

where R denotes a local matrix, then chiral symmetry is violated only by a 
local term 0(a). Dirac operators satisfying the GWC will be called GW Dirac 
operators. n 

LiischerEl has pointed out the explicit form of the associated symmetry of 
the action, 

ip -> cxp[i6»7 5 (l-a J RX>)] 0, -> ip exp [i (1 - a V R) 75] . (5) 

The GWC provides a sensible way (respecting the Nielsen-Ninomiya the- 
orem) to construct suitable Dirac operators. As has been shownEI such ac- 
tions cannot be ultra-local, but may be local, i.e. the coupling exponentially 
damped in real space. 
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fixed point action 



Figure 1. The FP action is the classical perfect action; the quantum perfect action may 
deviate from it away from the FP. 



For certain classes of Dirac operators one may derive stringent bounds 
on the shape of the eigenvalue spectrum. This, has been discussed in the 
framework of so-called fixed point (FP) actionsoO. 

On the lattice real space renormalization group (RG) transformations 
are realized through block spin transformations. The field variables over a 
localized region around a site x are averaged to produce the field variable 
of a coarser lattice. This mapping procedure in the space of configuration 
ensembles may be associated to a parameter flow in the space of actions. It 
is implicitly assumed that indeed a Gibbsian measure is suitable to describe 
the ensemble of blocked configurations. 

The continuum limit of the quantum theory is obtained at a critical FP 
of the system (cf. Figj^). The so-called renormalized trajectory in this space 
of actions leads to the FP and along its path there are no corrections to 
scaling. Ideally one could simulate the system at any point of the renormalized 
trajectory (with the so-called quantum perfect action) and obtain "perfect" 
continuum results. In reality that action may be quite complicated and has 
to be truncated in the number of coupling constants. 

Hasenfratz et al.El have suggested to determine instead the action at the 
FP of the RG transformation. For (asymptotically) free theories it may be 
determined from the classical field equations. The action has therefore been 
baptized "classical perfect action" . Its classical predictions agree with those of 
the continuum action independent of the coarseness of the lattice. FP actions 
are solutions of the GWC U. 

At the classical level, foii-FP actions, the Atiyah-Singer theorem finds 
correspondence on the latticeou; at the quantum level, no fine tuning, mixing 
and current renormalization occur, and a natural definition for an ordec, pa- 
rameter of the spontaneous breaking of the chiral symmetry is possible£3. R 
is then local and bounded and as a caascquence the spectrum of T> in complex 
space is confined between two circlescl (cf. Fig.||), 

|^ V-min | ^ Train : |^ Trnax \ — ^max > (6) 
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Figure 2. For FP Dirac operators the eigenvalues are confined in the shaded area of two 
circles tangential to the imaginary axis. The thick line indicates the position of the spectrum 
for R = 1/2. 



where the real numbers r m i n and r ma x are related to the maximum and min- 
imum eigenvalue of R respectively. For non-overlapping BSTs R = 1/2 and 
reduces to |A — 1| = 1, i.e. the spectrum lies on a unit circle. 
Independent implementations of the GWC are provided by the over- 
lap formalismEil, which allows the formulation of chiral fermions on the lat- 
tice. These, solutions are obtained in an elegant way, as shown recently by 
Neubergenlia, through some map of the Wilson operator with negative fermion 
mass. In this case we have R— 1/2 and circular spectrum |A — 1| = 1, too. 
For R = 1/2 the GWC assumes the simple form 

V + P f = D t £> = X>2> t . (7) 

On the lattice, an index for T> may be defined in a way analogous to the 
continuum, explicitly expressed by the relationa 

index(Z7) = -tr( 7 5 i?D) . (8) 

For GW Dirac operators and R—l/2 this relation comes out trivially. Only 
the modes with real non-vanishing eigenvalues contribute to the trace in the 
r.h.s; since the overall chirality must be zero, it reproduces up to a sign 
index(Z7). This index can be used to define a fermionic lattice topological 
charge, Qf e rm(U) = index (J7), for which the ASIT is satisfied by definition. 
For the FP action that fermionic definition—coincidesa with the pure-gauge 
quantity Qfp(C^), the FP topological chargeE3 of the configuration U: 

Q icm {U) = mdex(U) = Q Fp (U) . (9) 

The non-obviousness of this relation relies on the fact that Qfp(U) can be 
defined in the pure gauge theory, without any regard to the fermion part. 
This result is particular for a FP action and has no counterpart for a general 
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(non-FP) GW action. Of course, in practical implementations one relies on 
approximate parametrizations of the FP Dirac operator and the strictness of 
the relation is lost. 



1.3 Test Bed Schwinger Model 

As our test bed we consider a 2-dimensional (2D) quantum field theory with 
U(l) gauge group and Nf flavors of fermions. The action for the massless 
continuum model reads 



N 



S= I d 2 x [-F^F^+^fViPf] . (10) 



For Nf = 1 this is the Schwinger model0. 

This. 2D version of QED resembles 4-dimensional QCD in various 
waya 15 n 16 l. Quarks are "trapped" , i.e. in a mechanism superficially mimicking 
confinement we observe only bosonic asymptotic states. For Nf = 1 this is the 
Schwinger boson (called r\ by analogy to 4D) with the mass = g/n for the 
physical gauge coupling <?.._For the 2-flavor model one expects also a triplet 
of massless bosons (pions)LL3, although in 2D there can be not spontaneous 
symmetry breaking due to the Mermin-Wagner-Coleman theoremE3. Finally, 
there is a non-vanishing condensate (ip tf>) due to an anomaly. 

In the lattice formulation the gauge action in the compact Wilson- 
formulation is written 

S 9 = ^(l-Rc[/ P ), (11) 

where A denotes the lattice Zfj and the plaquette variable U p is the oriented 
product of links variables U x/1 G U(l) at site x in direction ^ = 1,2. In 
the continuum limit U x ^ ~ exp (i ag A x ^) where A is the gauge field in 
the non-compact continuum formulation, a denotes the lattice spacing, and 
/3=l/( 5 2 a 2 ). 

A lattice version for the integer geometric topological charge of the gauge 
field may be defined, 

Q(A) = -f [ d 2 xF 12 (x) Q(U) = -^Vlmlnf/p . (12) 

In torus geometryEl (i.e. periodic boundary conditions) this number may be 
non-zero for compact gauge field configurations. 
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The lattice fermion action is formally 



if) T) (m) tp 



(13) 



where T> (m) denotes the lattice Dirac operator matrix (fermion mass m) and 
the fermions are Grassmann fields. In a 2D context the Dirac matrices 7 M and 
75 are to be replaced by <r M and 03. 

In the quantized theory the Grassmann integration over the fermions 
yields factors (det T> ) and the expectation value of some operator C may 
be written 



where one usually samples over the gauge fields with some Monte Carlo pro- 
cedure. In the so-called quenched approximation (Nf = 0) the fermionic 
determinant is not included; this essentially neglects the fermionic vacuum 
loops. 

Dynamical fermions can be included either by incorporating them into 
the sampling probability measure for the gauge configurations or in the ob- 
servable. The first case is realized in the so-called "Hybrid Monte Carlo" 
method. However, since the determinant may be negative, only even num- 
bers of mass-degenerate fermions can be studied then. The second approach 
requires calculation of the determinant for each gauge configuration. It is 
plagued by the notorious sign problem, giving rise to possibly violent fluctu- 
ations and large statistical errors of the results. 

The gauge integral should sample over all topological sectors. For Dirac 
operators which allow exact zero modes the determinant weight (for massless 
fermions) removes the corresponding contribution to the integral. 

2 Lattice Fermions and Topology 

We have been studying the lattice Schwinger model for various lattice sizes 
and values of the gauge coupling, both quenched and with dynamical fermions. 
Here we discuss some of our results, with some emphasis on the properties of 
the spectra of GW operators. 

2.1 The Wilson Dirac Operator 

The original Wilson actionEl for fermions has the form 




(14) 



2?Wi 



m + 2) 1 M 

' 2 



(15) 
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Figure 3. Spectrum of the hopping matrix M for an 8 X 8 lattice at various values of /3, close 
to K c (f3); Q denotes the (geometric definition of the) topological charge. The eigenvalues 
for free fermions are degenerate. 



with the hopping matrix 

M xy = [( 1 + ^)U X yS x , y - ti + (l-a IM )Ut x S x , y+fl ] . (16) 

The hopping parameter n defines the fermion mass, for free lattice fermions 
K = 1/(2 m + 4). However, for the full theory with gauge interactions this 
relationship no longer holds. On one hand due to quark trapping there are no 
asymptotic fermions. On the other hand chiral symmetry is explicitly broken. 

Since this breaking is a local feature, we expect that at some K c {0) the 
chiral symmetry is formally restored in the sense, that the corresponding-Ward 
identity is satisfied with, a vanishing fermion mass parameter. This ideac3 has 
been pursued in QCDcj and we have utilized it also in the Schwinger model 
with Wilson fermions in order to identify the position of the critical line, where 
the fermion mass parameter vanishesE3. 

In Fig.^we compared typical eigenvalue spectra of the hopping matrix M 
at different values of the gauge coupling and in different topological sectors. 
We find strong correlation between the number of real eigenvalues and the 
topological charge. Closer inspection leads to the. conclusion, that indeed the 
real modes (counted according to their chiralitynj) may be interpreted as the 
would-be zero modes. Actually one has to divide that number by a trivial 
multiplicity factor of 4, since only the rightmost modes become zero modes 
of Pwi , whereas the other ones are doublers at other corners of the Brillouin 
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zone. The agreement of this number with the topological charge improves 
towards the continuum limit. 

In that sense even the Wilson Dirac operator obeys the ASIT for (3 — > 
oo; e.g. at j3 — 4 in 99% of the gauge configurations (sampled in a hybrid 
Monte Carlo simulation with dynamical fermions) the numbers of real modes 



(divided 
is obeyet 



4) agrees with the topological charge. Also the Vanishing theorem 



2.2 The Fixed Point Dirac Operator 

Although it was possible to obtain FP actions for some scalar or fermionic 
systems (cf. the summary in Ref.5) it turned out to be be a formidable 
problem for realistic gauge-fermion systems like QCD. We have succeedcdnJ 
to find explicitly an approximate FP action for the lattice Schwinger model. 
The FP Dirac operator is parametrized by a set of 429 terms bilinear in the 
fermion and anti-fcrmion fields, 

^Fp = £ Pi(f)*iU(x,f) (17) 

i,x,s 

The sum runs over paths / connecting a central site x with sites in a 7 x 7 
neighborhood; U(x, /) denotes the product of link gauge variables along that 
path. Standard symmetries are taken into account. Altogether we considered 
123 independent coupling parameters. 

The general idea is to start with a parametrization of the action, generate 
a set of gauge field configurations, perform the real space RG transformation 
and identify the blocked action. The matrix elements of the Dirac operator on 
the blocked configurations are then compared with the parametrization and 
the coupling constants adjusted. The whole procedure is iterated until the 
parameters converge. 

We studied samples of 50 gauge configurations and 14 x 14 lattices at large 
values of (3. We could show, that the resulting fermion action had couplings 
damped exponentially with their spatial extent; although the parametrized 
action is anyhow ultra-local by construction, this observation indicates that 
locality of an untruncated FP action seems feasible. 

Using the FP action we determined the bosonic bound state propagators 
for the Nf = 1 and 2 Schwinger model and found excellent rotational invari- 
ance and good scaling properties (cf . Fig J|) , to be discussed below in context 
of other actions. 

For the eigenvalue spectrum of the FP operator we found the situation 
shown in Fig.^: The eigenvalues are distributed close to a unit circle, in partic- 
ular towards the continuum limit (growing j3) . Towards the "hot" region the 
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Figure 4. Eigenvalues for sets of 25 configurations (superimposed) for lattice size 16 X 16, 
sampled according to the compact gauge action (from Ref.26). 



fuzziness increases. This demonstrates the predicted behavior for FP Dirac 



operators discussed in Sect. 1.2. The deviation from exact circularity is ex- 
plained by the approximateness of the action, which gives raise to fluctuations 
at smaller (3. 

Most remarkable is, that indeed we do find (within the numerical accu- 
racy) vanishing eigenvalues. These are individual zero modes and computing 
the matrix element (v^a^Vi) from their corresponding eigenvectors we find 
values ±1. Thus these zero modes have definite chirality. The partner modes 
with opposite chirality are at the other edge of the spectrum, at A = 2, irrel- 
evant in the continuum limit. The zero modes occur whenever the geometric 
topological charge of the gauge configuration is non-zero and its value agrees 
with the number of zero modes, even at finite (3. This agreement with the 
ASIT is quantitatively much better than for the Wilson action. 



2.3 Neuberger's Overlap Dirac Operator 

Motivated by the overlap actior0 Neuberger proposedlll to start with some 
Dirac operator with sufficiently negative mass, e.g. the Wilson operator at a 
value of m corresponding to ^ < k < 2 d-2 anc ^ then construct 

??Nc - 1 + 75 e(75 2?Wi ) ■ (18) 

Although the actual value of k used in this definition is largely arbitrary 
its choice may influence the approach to scaling in the continuum limit. The 
generalized sign function of the hermitian operator 75 Z?wi may be interpreted 
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Figure 5. The dispersion relations E(p) for the n (left) and the r\ (right) propagators, 
determined for three actions discussed in the text: Wilson (squares), fixed point (circles), 
Neuberger (diamonds) (from Ref.27). 



as a realization of 

l^wi I ~ ^/V Wi W Wi ~ 75 V(75^wi ) 2 

and is determined through the eigenvectors and the eigenvalues, 

e(75 V wi ) = U Sign(A) with 75 £> Wi =UAUK 

(Sign(A) denotes the diagonal matrix of signs of the eigenvalue matrix A.) 

Gauge configurations with non-zero topological charge implyo exact zero 
eigenvalues of 2?nc i therefore exact chiral modes are realizable. This was 
confirmed e.g. in Ref.27; we found an exactly circular eigenvalue spectrum 
and zero modes. Due to the freedom of choice of the hopping parameter in 
the Wilson action used to construct 2?Ne it cannot be excluded that some real 
doubler modes of the Wilson action are implicitly mistaken for zero modes 
leading to a violation of the ASIT at finite (3. This situation can be efficiently 
improved, if one starts already with a better action, e.g. an approximate FP 
action. 

An immediate question concerns locality of this operator. In 4D QCD 
it has been demonstrated!^ that for large enough (3 one may expect locality 
and for a Monte Carlo generated ensemble of gauge configurations exponential 
falling off of the effective coupling parameters has been shown. 

Neuberger's operator provides an explicit example of a GW action with 
optimal chirality properties. However, nothing can be concluded as to the 
scaling properties towards the continuum limit. In fact, in Fig. ^|we see, that 
the spectral properties of the bound states propagators are not improved over 



(19) 
(20) 
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those of the original Wilson action. Pno is automatically 0(a) corrected!!!! 
and thus improves scaling for the on-shell quantities; without at least introduc- 
ing improvement of the current operators one would not expect improvement 
for the propagators, as shown by our results. The two features - chirality and 
scaling - seemingly may be considered quite separate issues. Optimal actions 
should allow chiral fermions and have good scaling properties. 

2.4 Spectral Distribution and chRMT 

The limiting value of the spectral density for small eigenvalues and large 
volume, 

- 7T lim lim p(X) = (ipij>) (21) 

provides_^an estimate for the chiral condensate due to the Banks-Casher 
relational. Exact zero modes are disregarded, only the density close to zero 
is of relevance. 

In order to study this property of the lattice spectra (for the GW Dirac 
operators) one should map the eigenvalues from the circular shape to the 
imaginary axis. This is done by a stereographic projection 

^ = A- (22) 
1 2 

Near zero the resulting distribution on the (tangential) imaginary axis agrees 
with that on the circle. The generalization of such a projection for a general 
GW operator is the introduction of 

t = Y^RV ™ 
which anti-commutes with 75, is anti-hcrmitian and has purely imaginary 
spectrum. A procedureJo define a subtracted fermion condensate has been 
suggested by Hasenfratza (within the overlap formalism cf. Ref.32), 

\ \ / / gauge \ \ *s / / gauge 

(see also Ref.s 12,26,32,33). The expectation value includes the weight due to 
the fcrmionic determinant. 

In Fig.^| (from Ref.27) the spectral density for the Neuberger operator is 
compared with that of the FP operator; near zero we find good agreement 
with each other and with the value expected from the continuum theory. 

Studying the spectra, of the Dirac operators suggests comparison with 
Random Matrix TheoryE£l (RMT) . The spectrum is separated in a fluctuation 
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Figure 6. The unquenched (Nf = 1) eigenvalue density distribution (16 X 16, (3 = 6) 
projected from the circle onto the imaginary axis for fpp (thick lines) and X>N e (thin 
lines). The horizontal line denotes the continuum value at infinite volume. 



part and a smooth background. The fluctuation part is conjectured to follow 
predictions ixL.one of three universality classes. For chiral Dirac operators 
(chiral RMTE3) these are denoted by chUE, chOE and chSE (chiral unitary, 
orthogonal or symplectic ensemble, respectively). Several observables have 
been studied in this theoretical context. Comparison of the data should ver- 
ify the conjecture and allows one to determine the chiral condensate. This 
information is contained in both, the smooth average of the spectral distribu- 
tion, and in the fluctuating part. In particular the distribution for the smallest 
eigenvalue p mm (A) contains this observable: Its scaling properties with V are 
given by unique functions of a scaling variable z = XVE, depending on the 
corresponding universality class. Usually this is the most reliable approach to 
determine S, which then serves as an estimate for the infinite volume value 
of the condensate in the chiral limit. _ . 

Within the Schwinger model we founcE-J the universal properties of the 
(expected) chUE-class, unless the physical lattice volume is too small. In 
Fig. ^ we show the distribution of the smallest eigenvalue for two of the 
actions studied. Further examples, al^-iaxJiiding dynamical fermions and 
4D applications, are discussed in Ref.soEZlES. It has become evident, that 
chRMT describes the fluctuating part of the spectrum satisfactorily within the 
expected universality classes. It provides a means to separate the universal 
features from quantities like the chiral condensate, which have dynamic origin. 
Or, to say it (tongue in cheek) more provocative, to separate the known thus 
uninteresting, universal properties from the unknown, physical properties of 
the system. 
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Figure 7. The distribution density for the smallest eigenvalue in the quenched (Nf = 0) 
case in the sector with topological charge zero. The full line gives the prediction for the 
chUE class, the broken line would correspond to the chSE universality class. 

3 Conclusion 

The main messages I wanted to deliver here are 

• Lattice gauge theory has made a big step forward in the understanding 
of chiral fermions and their lattice formulation. The Ginsparg- Wilson 
relation is a cornerstone in that development. Meanwhile U(l) chiral 
gauge theories with anomaly-free multiplets of Weyl fermions have been 
constructed properly^. 

• Neuberger's overlap action is a wonderful testing ground for chirality; its 
implementation in 4D is computationally expensive. Perfect actions ide- 
ally provide the best of both, chirality and scaling, but their construction 
in realistic 4D models seems to have unsurmountable problems. In a 2D 
gauge theory, however, we were able to construct an explicit example, 
which has all the beautiful features expected. 

• The spectra of Dirac operators have universal features, which are de- 
scribed by chiral RMT. This opens the path to efficient methods to ex- 
tract from the data non-perturbative quantities like the chiral condensate. 

• The Schwinger model has provided an excellent test case to study various 
features also observed in QCD. Although maybe of different dynamical 
origin, mechanisms like confinement, chirality and symmetry breaking 
can be nicely studied in their lattice realization. This may lead to a 
better understanding also for the 4D theories. 
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